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On Nearly Triangular Matrices 



A. M. Ostrowski 2 

A discussion is presented of the change in the inverse of a triangular matrix if on one 
side the zeros are replaced by sufficiently small numbers and on the other side the non- 
vanishing elements are varied by sufficiently small amounts. 



1. Introduction 

Consider a system of linear equations 

n 



(1) 



with the matrix A, in which all diagonal elements a^ (m=1, . . . f n) are not equal to and 
the elements off the diagonal satisfy for two positive numbers m, M the inequalities 



VI = m \ a n, 



K"l ^M\a^ 



(j><ji;M=l, • . ., n), 
(?>/*; m=1, . . ., n-\). 



(2) 



If m is very small, the system does not essentially differ from the corresponding "triangular" 
system in which all a M „ with *><Cm are replaced by zeros and the matrix of which will be denoted 
by A m . It then appears plausible that the solution of this triangular system does not differ 
very much from that of the system (1). 

However, the value of the determinant of the order n 



-M 

1 




-M 













-m 







1 -M 

1 



= l—mM n -\ 



shows that if M is, for instance, greater than or equal to 10, the determinant of our system 
will not be necessarily different from zero unless m< 10" (re_1) . A detailed study of the problems 
connected with the matrices characterized by (2) appears, therefore, to be of importance and 
interest. 

As the first problem in this connection, we give a necessary and sufficient condition that any 
matrix A satisfying conditions (2) be nonsingular. If m<^M, this condition is given by 



in 



(l + m) n ^-(l + M) n ' 



(3) 



i This paper was prepared under a National Bureau of Standards contract with American University. 
2 American University and University of Basle, Switzerland. 
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and, if m=M, by ; 



m< -• (3°) 

71—1 X ' 



In order to obtain a precise measure of the influence of the change from A to A {Q \ we have to 
discuss the estimates (for convenient norms) of the norm of the matrix A~ l — A (0) ~ l . 

We consider, in particular, two such norms defined in section 5 and denoted by 
\A~ X — A^^lp, (p=l,co), which are particularly suitable for the problems of numerical analysis. 
Assuming, without loss of generality, that a tlfl =l (m=1, • • . , w), we show that for given 
values of m, M such that M2£ 1.5/Vfc, n^4, we have 

\A- l -A™- l \ v ^(l + M) n - l -^~ (M^—,n^4), (4) 

1 — a \ n / 

where 1 — 6 is the smallest modulus of the determinant attainable for the matrices A and is 
connected with m by the relation 

5 =J^L, <«1, (5) 



with 



TO 



M ^ (1 + M)«-.M-1 , (6) 



If M<0-5/w, the formula (4) need not be valid any longer, but we can prove in this case 
the relation 

]A -i_ Am - ilp ^_^m_ (M£1.5fn) (7) 

is valid as long as m remains less than l/2n. 

The estimate (7) is not a "best" estimate for all values of M^1.5/n, but still it is not far 
from the best, since for m=M^l/(n—l) we have 

l^-i-^(o)-i| < fr-i)"* ( m =M<~--\ (8) 

lp 1 — (n— l)m \ Ti—1/ 

which cannot be improved for any value of m<Cl/(n— 1). 

The condition (3) is derived in section 4, theorem B. However, we derive it as a special 
case of a more general theorem, where in the inequalities (2) the expressions m, M depend on /x, 
that is to say, change from one row to another. The necessary and sufficient condition for all 
matrices A to be regular (theorem A, section 4) is in this case rather unwieldy, but still may be 
very useful in some cases because it contains 2n— 2 instead of two essential parameters. The 
direct derivation of theorem B is, of course, much simpler, since, as the reader will immediately 
see, the computations of the determinant n in section 2 can be considerably shortened in this 
case. The connection between the formal algebra of sections 2 and 3 and theorems A and B 
is provided by a result concerning the so-called ^-determinants and M-determinants published 
16 years ago [1]. The results about the norms I^L" 1 — A m ~ l \ p are obtained by using the 
explicit representation of the inverse matrix of a certain matrix A w , which provides a majorant 
for all matrices A~ l . The formulas giving A" 1 are derived in the second part of section 3, and 
in section 6 the norms |A~ ! — A^" 1 ^ are derived and discussed. The corresponding inequalities 
for |A _1 — ^l^" 1 ^ are then obtained in section 8, by using a new theorem (lemma III) concern- 
ing the connection between the i7-determinants and the M-determinants. 



3 This condition (3°) is already contained in some results in a paper [2] 4 , and also can be deduced from a well-known theorem of the theory of 
determinants discussed in [4]. 

4 Figures in brackets indicate the literature references at the end of this paper. 
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In this section is given another application of this theorem in estimating the variation in 
the inverse matrix of a triangular matrix satisfying the conditions (2) with rn = 0. We obtain 
an unexpectedly simple and elegant formula (120). 

In section 8, the results are applied explicitly to the problems concerning* the linear system 
(1). It may be finally remarked that these results, with obvious changes, remain valid if in 
the matrix of (1) the rows and the columns are interchanged, although no mention is made of 
it explicitly at every step. 5 



2. Value of the Determinant 12 n 



Let K n be defined by 



K n 



K n 8 n 

K n _i K w _i 5„_! 
K n __2 Kn-2 K-n-2 &n-2 



K2 



Kl 



*2 



Ki=k u K 2 = 



*2 

K 2 5 2 










K 2 5 2 
Ki Kl 

= ki(k 2 — 8 2 ), 



(71^3), 



(9) 



where in the /u-th row all elements to the left of the main diagonal and on this diagonal are 
equal to /c w _ M+ i, the next element to the right is 5 w _ M+ i, and all other elements are 0. k m , 5 m 
are here independent variables. In subtracting the second column from the first, we obtain 
K n =(K n —8 n )K n _i, and therefore the following formula, \alid also for n=l, 2: 



K n =Ki n (k v — 8 v ). 

v=2 



(10) 



Consider now for n ^ 3 the determinant 

b n 
7 tt _i 5 W _! 

7 n _2 Ytt-2 $n-2 

T = 



72 72 72 72 

1111 



8 Some of the results contained in the sections 2 to 5 have been published w 



o p, 
o p„-i 
ft,., 



s 2 ft 
1 

ithout proof in [3], 



(»£3). 



(ID 
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We have in particular 



T,= 



B s ft 
72 S 2 ft 
1 1 



■(5 3 ft + ft(52-72)). 



(12) 



Developing T n in the elements of the first line and using the value (10) of K n , we obtain 
for n > 4 



T n = 8 n T n -i — f3 n IT (5,-7,), 



T 



Tn-l ft, n Z 



&2 • ■ ■ S n 62 . . . 5„_i S 
and therefore generally for n ^ 4 



„=2 \ 8 V J 



8 2 . . • 8 n 

Since by (12) 



we obtain 



P'nYi ^i^-'nYi 7> ")| ... + ^n (i_2*Y] + A. (13) 

• • • 8 n M = 2 5 M „ = 2 \ £„/ 



A. 

^2^3 



where II is identically 1, and therefore finally 

*=2 



Tn^-itp* n 5, n 1 (5,-7,). 

M =2 v=/i-(-l v=2 



If we now put 



T* = 



8 X 

72 5 2 



Tn-l Tw— 1 Tw— 1 

1 1 1 



o ft 
o ft 



Tn-l fti-i 
1 



(n^3), 



21= 



5, ft 
72 5 2 ft 
1 1 



(14) 



(15) 



this becomes T n if the indices of /?„ y v , 8 V are replaced by their complements with respect to 
n+1. We obtain then from (14) 

n n n—1 

t*= — S ^n+i-M n ^w+i-v n (6 w+ i_ v — 7 n+ i_ y ), 

M=2 y=/x+l v=2 
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or in replacing the summation index jjl by n+1— k 



Tn =—^2i Pk n 5 n+ i_„ n (5 n+ i_ v — 7 W+ i_ v ) 

k=1 v=n+2—K v=2 

and finally, if in both products v is replaced by n-\-l — X, 

W — 1 K — 1 71 — 1 

T^-S&nsx n (5x-tx). 

K = l \ = i X = /c+l 

Consider now for n^3 the determinant 

ai —Mi —Mi 

—m 2 a 2 —M 2 

— m z — m 3 a 3 



I", 



-m n 



-m„ 



-M 2 
-M 3 



If we subtract here the last column from each of the preceding ones, we obtain 
Mi + ai ... —M x 



M 2 —m 2 
M 3 —m 3 



A'f 2 -{-a 2 
M 3 —m 3 





M 3 + a 3 



M n _i—m n „i M n _ t — rrin-i M n _ x — m n _i 
-(a n +m n ) -(a n +m n ) —(a n +m n ) 








-M 2 
-M 8 



M n -i— m n _i -M n _i 
— (a n +m n ) a n 



(16) 



(17) 



Here the subdeterminant corresponding to the last element of the last row is obviously 

n—l n—1 

n (My + a,), so that ft w is the sum of a n IT (M„+a„), and 



(a n +m n ) 



Mi+ai 








M 2 —m 2 


M 2 -f- a 2 





M 3 —m 3 


M 3 —m 3 


M 3 + a 3 



M n _i—m n _i M n _i—m n _i M n _i—m n _i 



1 



1 



1 






-M, 





-M, 





-M, 



1 
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This last determinant becomes T n * if we put 

b v =M v +a Pj y v =M P —m v , $ V = — M V (v=l, . . ., n— 1), 
and has therefore by (16) the value 

— S(-M«) n(M x +<*x) n (wx+ax). 

K = l X = l X = /C+l 

We obtain therefore for ft n the expression 

71 — 1 71 — 1 K — 1 71 

Q>n=<xn n (M„ + a„) — SM» E (Mx+ax) n (^x+ax). 

V = l K = l X=l X = K+1 

The determinant tt n can be also written in the form 

ocn —m n —m n 

—M n _ x a n ^ — m n _! 
-Mn-2 — Af n _2 a w _ 2 



fi„ 



-M 2 — M 2 —Mo 
-Af! -M! -M x 



—m n —m n 

— m n _ x — w w _! 

— W n _ 2 — W w _ 2 



a 2 —^2 

-Mx a, 



(18) 



(19) 



and we obtain therefore from (18) 



7i n 



Sl n = aiU (m v +a v ) — X>* II (™m+«m) n (M M +a„). 

^ = 2 K = 2 m = k+1 M=l 

3. The Matrix A n and Its Inverse 
We consider now the matrix 



r 1 -M -M 

-m 1 — M 



AT| 
-M 



. A n = 



v.— m — 7/z — m ... 1 ./ 
Its determinant is obtained from fi tt in putting in (20) 

ai = ot2= . . . a»=l; m 2 =rn 3 = . . . =w n =m; M X =M 2 = 
We obtain 



(20) 



(21) 



=M JI _ 1 =M. 



324 



and this becomes, if M?±m, 



| i .|-(™ + l)--».(M-H) CM+D^ + l)-- 

[(m + l) w - 1 (M-m + m(M+l))-m(M+l) n ], 



~M—m 



\*n\=j^[M(m + iy--m(M+iy] (m^M), 6 (22) 

while for m=M we obtain from (22) in letting Al—>m: 

|A n | = [l — (tt--l)m](l+w) n "" 1 (rn=M). (23) 

In particular, if 0<ra<M, a necessary and sufficient condition for |A n |>0 is 

m ^ M , . 



(ra + l) n ^(M+l) w 
We assume now in particular 

0<Cm<M. (25) 

If we introduce the abbreviations 

(l+M)*-nM-l (l+ro)»-nm-l 

Mn= w ~' Mn== ~ ~m ' C25) 

we can write (22) in the form 

(M—m)\A n \=M(l+nm+mm n )—m(l+nM+MM n )=M—m—mM(M 
and therefore, if we put 



d=l-\A n \, 






(27) 


8 M n — m n 
mM M — m 






(28) 


It follows from (25) for n>2 








M n -m n * (n\ M v - l -m*- 1 ^ ^ (n\ , 2 
Af-ro £*W M-m ^&W 


(1 + M)*— 1- 

M 2 


-nM M n 
~ M' 




so that from (28) we have 








M " > 5 >M n . 
mm 


- 




(29) 


l ~M 









It follows in particular that if (24) holds, then 5>0, 0<|A n |<l. 
In solving (29) with respect to m and d, we obtain 



M„m< 5 <^f, l-J^, 0<|.<1. (30°) 



6 This formula can be also obtained from the formulas given in the proof of theorem III in [3a, p. 1131. 
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In order to discuss the meaning of the condition (24), consider the curve 



we have 



(l+x) n 
,,, , \—{n—l)x „, N (n— l)x — 2 



(31) 



(l + a)" 



(l+^) w+2 



For a; >0 the curve (31) has one maximum at ic= 1/(71—1) and an inflection point at a; =2/(n— 1). 
In figure 1 the curve (31) is drawn (computed for n=5). The portion of this curve from 
the point to the highest point T will be denoted as the ascending branch and the portion 
between T and x=°° as the descending branch. 

f(x) 



I I 
j I 



m m(M) -±- 



M JL 
n-i 



f(x) = 



(i + x) 5 
Figure 1. 

If we assume that in (24) m is less than M, we see that either ■ 



0<m< 



n—1 



My 



n—l 



or 



m<M^- 



(32) 



(32°) 



In order to find for a given M^>l/(n— 1) the range of values of m satisfying (24) we find 
on the curve the second point P M with the same ordinate as the point x=M; then if the abcissa 
m(M) corresponds to P M , we have (X^nK^miM). 

Finally we prove from (24) and m<CM that it always follows that 



mM n ' x <y. 



(33) 



Indeed, in proving (33) we can obviously assume that M is greater than 1 and therefore 
on the descending branch in the diagram. But then we have 



1 



M n - 



1 + 



1 1 + M'> 



w-l /" 



M 

l + M' 



M- 1 

and therefore in raising it to the nth power 

1 



M»- 



( 1 + i^) 



»> 



M 



(1 + M) n 



We see that l/M"' 1 must lie, in figure 1 , between M and m(M), and it follows that l/M B_1 >m, 
that is (33). 
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We are now going to determine the inverse matrix 



(34) 



of A», if its determinant (22) resp. (23) is not equal to 0. It is sufficient for this purpose to solve 
the corresponding linear system 



M — 1 7? 

Xft—m^Xy—M Xj x v =z it (jjt=l 9 . . .,n) 

V = l V = fl-\-l 



for indeterminate z„. Put 



Af 



m — M 



1 + AT 



'X, 



1+m 

1 + (r= l + M = ^ 



a,, + \s = ?/„ 



Then the system (35) is equivalent with 

(1+M)x fi =(m-M)(x l + . . . +x^ l )+Ms+z ll (ju=l, . . .,n), 
aj M =<r(a?i+ . . . +.r M _,)+^ (m=1, . . .,^), 

S v =x { +. . . +x„ (p=1, . . . ,ft), S =0: 

Sn—Sp-^eSt-i+yM S„=(l+a)S^i + y ll =qS l ^i + y ll9 

S.q-^S^.q-^ + y.q-^ (/i=l,. . .,»), 



or in putting 



(35) 



(36) 



S>=S 2/"^"" (M=l,*. • .,n), 



M-l 



*M = S^(?-l)g' l "'- 1 + »M (M=l,. - .,*). 



From (37) for'/i=n we have by (36) 



J>=1 f= 1 v= 1 



q— l v=i 
where for m=M, q=l the coefficient of \s is to be replaced by n. Since for m t^M 



q n —l _ m — Mq n _ Mq n — m 
q — 1 "" m — M ' M—m 



we have from (39) and (22) 



s=QSM w " y (m^ili), 



Q* 



M—m _(M+l) n 
~Mq n — m~ |A B | 



(m,9*M), 



(37) 
(38) 



(39) 



(40) 
(41) 



294330—54 5 
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while for m=M we obtain again (40) in defining Q by 



n= m + 1 - 

^ 1— (n— l)ra 



|A TO | 



(ra = M). 



If we now replace in (38) y v by \s+a v and use for s the expression from (40), 
av=(ff-l) S a^-^+^+Xs P(g— 1) S ff'-^+lT 
But here the expression within the brackets is <z M-1 , and we obtain therefore 

and finally in introducing a M from (36) 

^-§^'" , ^+3fti+^$rli «'«■-; 0.-1,... ,»). 

Introduce the two following triangular matrices: 

ro i 



U n (q) = (u lxv ) = 



1 . 
q 1 
g 2 q . 



<z 2 <z i oj, 



GO") 
(mS*), 



f0 g- 1 q~ 2 q~ z 



g" 1 q~ 2 



~-w+r\ 



^»fa) = (v) = 



J) ... 
We see then that in (42) the matrix corresponding to the first right-hand sum is 



V- 



q-1 _ m-M 



M+1 "" (M+1) 2 



(41°) 



(42) 



(43) 



(44) 



(45) 



(46) 



328 



whereas to the second right-hand sum in (42) corresponds the matrix 

where E n is the unit matrix. We obtain 

/ m-M MQq" \ ( 1 > MQq«^ \ MQq- 1 

A " V(M+l) 2t (M+l) 2 r wt U+l t (M+l)V Bt (M+l) 2 Vn ' 

The coefficient of U n is here equal to (mQ)/(M+l) 2 , both in the case (41) and (41°), and we 
have finally 



A -i = m Q 

n {M+iy 



(7n+ U+l + (M+l)V w + (M+l) 2 Vn ' ^° 



We now denote by s M the sum of the elements in the /xth row of the matrix (47) and by £ M 
the corresponding sum in the /xth column. If we first assume m^M, for the matrices U n 
and V n , the fith row sums are, respectively, equal to 

g"" 1 — 1. 1 — q ti ~ n 

1 — ' 1 — ; 

q-1 q-1 

and the *>th column sums are 

q n - v -l . , l-g 1 - y 
q-1 ' q-1 ' 

We obtain now from (47) by (36) 

(M+l) 2 s,= mQ^ l ^ + M+l+MQq n - 1 + MQ qn ~ 1 ~f~ 1 
q— 1 q— 1 

=^[mq^- 1 -m + Mq n -Mq n - 1 + Mq n - 1 -Mq^ l ] + M+l 
cr 

1 Q [Mq n -m + (m-M)q fi - 1 ] + M+l f 



m—M 

and this is by (41) equal to (M-f l)(?g M_1 . We have finally 



^X r iJM+ltl H (48) 

^~M+1 3 |A n | 2 ' l48j 



and this remains true also for m=M } as is immediately seen. 
Similarly, we have 



n -1 ^n -v 



and we have 



(M+l^-^^^^M+l+M^g-^ + MQ ^"^ 

0" 

t= ^M Q[M9 " _m+(m_M)2 "" r]+M+1=(M+1)Qr "'' 



** Af+1 2 |A.| 9 ' [ } 



relation which is also immediately verified for the case m=M. 
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4. Bounds for the Determinants, depending on &n and A n 
A determinant 



M= 



Oil 

-m 21 



-m 12 
a 2 



-m rt 



-m n2 



— mm 

— m 2n 



(50) 



will be called an M-determinant if all diagonal elements a v are positive, all elements off the 
main diagonal —m^ (ixt^v) are not positive and the determinant M, as well as all principal 
(coaxial) minors of M of all degrees, are positive. 

In what follows we will have to use a theorem given by the author in 1937 [1], and which 
will be formulated as the following: 

Lemma I. If we have for the M-determinant (50) and a determinant H=\h ltv \ (n,v=l, . . ., n) 
of order n, the inequalities 



\h vv \^a Vi Ih^l^m^ (n?*i>; jjl,v = 1, . . ., n), 



then H^O, and we have 



\H\^M. 



(51) 



(52) 



If (H), (M), respectively, denote the matrices of the determinants H and M, the inverse matrix 
of (H) is majorized by the inverse matrix of (M) 



{H)- l <.(M)-\ 



(52°) 



Suppose now that in the determinant tt n given by (17) all a v are positive and m v , M v non- 
negative. Then it follows from (18) that Q, n is a monotonically decreasing function of all m v 
and from (20) that M is also monotonically decreasing in all M v . Suppose now that for a 
certain set of values of <x v , m v , M v the determinant £2 W ^0. Then replace the m v and 
M y corresponding to certain rows (v ly . . .,v K ) by zeros; the determinant Q n cannot decrease 
and remains therefore positive; but then ti n becomes equal to the product of a Vl a l>2 . . . a v< 
with the principal minor complementary to the set of indices v u . . .,v K . Therefore, all principal 
minors of 12 w are positive. We see that 12 n is an M-determinant , if the conditions 



a v >0, m v ^0, M„^0 (*=1, 



,n); 



(53) 



are satisfied and 12 w ^0. But now we can easily deduce the following theorem. 
A. Consider the set of all determinants ^4=|a M „| satisfying the conditions 



M^m^ (*</*), M^M* (">m), 



I^MMl = a M 



0v=l, . . . ,ft), (54) 



where a M are n given positive constants and m M , M^ are 2n—2 given nonnegative constants. Then 
in order that no determinant A of this set vanishes, it is necessary and sufficient that ft n >0. If 
this condition is satisfied, we have 

Proof: Since in the case 12 M >0, 12 n is an M-determinant, the sufficiency of our condition 
follows immediately from the inequality (52) mentioned above. 



If tt n =0, we can take obviously a M 



-dp and a^- 



-m u 



or a yLV =—M tl according as 



*></z or j>>jLi and obtain a vanishing determinant A satisfying the condition (54). 
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Suppose now £2 n <0, then it follows from the form (17) of tt n that ft n becomes positive if 
all m M are replaced by zeros. There exists, therefore, such a positive £<1 that Q n vanishes if 
all m M in (17) are replaced by tm^, but then we obtain a vanishing determinant A of our set 
in taking a llll =a lt and a flv =—tm fi or a fJLV =—M^, according as i></z or *>>/*. The theorem A 
is proved . 

In specializing the matrix of il n to A n and in assuming that in particular 0</ft<M, we 
obtain immediately from (22) and the theorem A: 

B. Consider for two positive constants m, M, m^M, the set of all determinants A=\a fiy \ ojthe 
nth order for which 

|aj£l (m=1, • • • ,»), M^w (*</*), k.,|^M (y>/i), (56) 

fAm, m order fA#£ no determinant A of this set vanish, it is necessary and sufficient that for 
m<CM the inequality 

(ifb*<(iw» {m<Mh (57) 

holds j and if this inequality is satisfied, we have for each determinant A of the set 

Mima < q+^) w (i+Mr r M m 1 

1-AlfciA.i- M _ m l(i + M) n (l+ro)"J [bb) 

If m = M, the condition becomes 

m<r~zzi (M=m) (57°) 

and we have, if (57°) is satisfied, 

\A\^\A n \ = [\-(n-l)m](\-\-m) n - 1 (m=M). (58°) 

From theorem B we can deduce the following theorem. 

C. Let A={a flv ) (n,v=l, . . . ,n) be a matrix satisfying the conditions 

\a„\&m (*<», M^M (y>ii), (59) 

where m, M are two constants with 0<ra<M. Put 

8(m,M)= M Mn_ m Mn ' (60) 

7%eft all fundamental roots of the matrix A are contained in the set of the n closed circles described 
around the elements a^ with the radius 8(m,M). The value (60) of 8(m,M) cannot be improved 
ifA=A n . 7 

Proof: Let X be a fundamental root of A so that the matrix \E—A is singular. Put 
min|X— a^\ = a; we then have to prove that a^8(m,M). If a=0, there is nothing to prove. 

A* 

Suppose a>0, and consider the matrix 

For this matrix we have 

|*Jfcl, IU^ (Km), \K,\^ (?>*). 

a a 



1 For M=m f 8(m,M) becomes (n— \)m, but in this case our result is contained in Gerschgorin's theorem, cf. 14J. 
From the theorem C the result given by Stein and Rosenberg in L3aJ as theorem III follows immediately. 
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Therefore, we have by theorem B 

m M 



a a 



and therefore, 



i+= i+M 

a a 

~wF n ~ =z ~W TrrJ 

( 1 1 \ Af m 

°\™ 1/n M 1/n / = M 1/w /n 1/n ' 



^ m -Un_ M -lln -Km,M). 

5. Bounds of the Matrix A^ 1 

For an n dimensional vector £=(#!, . . . ,x n ) the Holder norm corresponding to the 
exponent p ^ 1 is given by 

|{|, = VI*l| f + • • • +\*n\ P (PSD- (61) 

We will only use the three cases corresponding to p=l, 2, oo : 



|€|,='W+ • • - + W> |*L=max|s,|, |£| 2 =VW 2 + . . . + W 2 - (62) 
We have among these three norms the following inequalities: 

|{|.S|*|i£»|*L, (63,a) 

|«I.^l*|f^V»l€U, (63,b) 

|«|iSVn|«|iSV»l«k, (63,c) 

which are immediately verified. The left-hand inequality (63, c) implies the well-known 
inequality between the arithmetical mean and the arithmetical mean of the squares. If 
A= (a M „) is an n Y^n matrix, we define its norm corresponding to the exponent p (1 ^p g oo ) by 

\A\,=max^ (64) 

and denote it by \A\ P . We will use here too only the cases p= 1, 2, oo . 

In applying to the definition (64) the formulas (63, b), and (63, c) we obtain immediately 

-L\A\ t £\A\,*Jn\A\u (65,a) 

±\A\.£\A\ a £,/n\A\„. (65, b) 

For p= 1, oo the expressions of \A\i, \A\ m are easy to write down; we have, as is well known and 
very easy to prove, 

|^4.|i=maxX/ \ a n>>\ (66, a) 

|il| 1B =maxS W- ( 66 »b) 
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As to \A\ 2 , its expression is irrational; \A\o is the square root of the maximum fundamental 
root of the symmetric and nonnegative matrix AA*. Since the direct computation of \A\ 2 is 
difficult in most cases, we prove the following estimate for \A\ 2 : 
Lemma II. We have for any matrix A 

imaxd^W^IJ^I^I^VI^Iil^l., (67) 

lb 

The first part of (67) follows from (65, a) and (65,b). To prove the second part, we introduce 
the notations 

s M =X) \ a ^\y iv = ^2 |<v| (/i,p=l,. . ,,n). 

v n 

The sum of the moduli of all elements in the /zth row of A A* can be estimated as follows 

V, K K V K 

Equation (67) follows now from the theorem of Frobenius that the modulus of each funda- 
mental root of a square matrix does not exceed the greatest sum of the moduli of the elements 
of this matrix in different rows. 

If we now apply these results to the matrix A" 1 discussed in the section 3, we obtain from 
(48), (49), (66,a) and (66,b) 



and from (67) and (65, a) 



|A„ |.-|A„ U-jyqrj- — ^ — » (68) 

Q <| A -»i < Q _ (^+D"- 1 (m 



In combining these inequalities with the results given in section 4, we obtain the following 
theorem. 

D. Let A—ia^y) be a square matrix of order n satisfying conditions {56), and let (57) be 
satisfied. Then we have 



|A-|,.^±fl' (,= 1,2,.); 



and therefore for any vector £ 



l^*l»^ (M+i 1 )-' 1 * 1 ' (P =1 > 2 > ro )- ( 70 > 

Proof: It follows from (52°) of lemma I at once in virtue of (66,a) and (66,b) that 

lA-^^A^^j^- (p=l,.«0; 

further, the matrix ^L -1 ^ -1 )* is majorized by A" 1 ^" 1 )*. By a well-known theorem of 
Probenius, therefore, the maximum modulus of a fundamental root of A~ l {A~ 1 )* is majorized 
by that of A" 1 (A" 1 )*, and so we also have 

By definition (64), therefore, in putting Ai£=ri. w e have 

and this is equivalent to (70). 
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6. Bounds of the Matrix A^-A* ' -1 
We denote by A™ the matrix obtained from (21) in replacing m by zero: 

r l -M -M . . . -M' 



A C0) = 



(i 



-M 











-M 



1 J 



(71) 



We will prove the following theorem: 

E. If 0<m<M and (57) is satisfied, then the matrix A" 1 — A^" 1 is a nonnegative matrix 
and we have 

|A.-'-A«°>-*l,,=max (1 + M)- 1 ( 1 + W )'~'-I A »I (p=l,co). (72) 

Proof: Since A™ satisfies the inequalities (51), it follows from (52°) that A" 1 — A"" _1 is 
nonnegative. Denote by s^ and ^ 0) the sums of the elements in the /xth row and in the *>th 
column of A®*" 1 , and by s M and t v the corresponding expressions for the matrix A" 1 — A^" 1 . 
We have 



VfS/i— $£ 0) , t„=t v —ti 0) Gu,p=l, . . . , n). 

It follows then from the formulas (48) and (49) that 5 M runs for fjL=n, n—1, 
same set of values as t 9 for v=l, . . ., n. 

Formula (49) can be written in the form 



fcn 



g w ~ y (M+ l) n - x = (l +m) w - y (l + M) y ~ 1 



(73) 
1 through the 

(74) 



Since A„ becomes 1 for m = Qj we have (72) and E is proved. 

Discussion of [A" 1 — A^" 1 ^. We prove first that t v — t, i0) increases with v as long as the 
condition 



i-S^i 



is satisfied. Indeed, put 



k>= T /i"' =(M+l)- ( ' 






(l+m)»--|A n | 



In solving the three inequalities 



k,= l 



(75) 



(76) 



(77) 



with respect to |A n |, we obtain correspondingly 

| As || ff , s ( 1 _^(l +m )«— i. (78) 

The inequalities &„ ^ 1 (v—1, . . ., n) are obviously satisfied in virtue of (75) as long as 
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v^n—1. We have, therefore, in this case, in using (72) 

Max {t-in = (l + M)- 1 y^W + M)"- 1 j^-, 

v |A n | 1—0 

|A- 1 -Ar- 1 | 1 , = (l+M)"- 1 ^^=(l+M)»- 1 r ^ (p=l,°>; [Anl^l-^y (79) 

Condition (75) can be written by using the notation (27), as 5^m/M, and in virtue of (29) this is 
certainly satisfied if we have MM n ^. 1, >• 

(Al+l) n -nM-2^0. (80) 

Condition (80) becomes for n=3: M^O.5321. For ti^4, (80) is in any case satisfied, ij we have 

M>— (ti>4). (81) 

~~ n ~ 

Indeed, if the relation (80) is satisfied for a positive M, it is satisfied for any greater value 
because the coefficients of all positive powers of M in the left-hand expression are greater than 
or equal to 0. To prove the sufficiency of (81), it is sufficient to prove that (1 + 1. 5/71)* 2s 3.5 
(n^4). But here the left-hand expression is mono tonic ally increasing with n, and this 
inequality follows, therefore, from (5.5/4) 4 =3.75 . . . >3.5. 

Suppose now that we have 

1-5*1*1*1. »*£■ (82) 

Then we have, since 0<m<l/(7i— 1), 

s±^ i K^ s( ,_ 1) (( 1+; i I j-'_ 1 ) s(1 _ 1)( ,, ) , 

71-1 

(l+m) w - y =l + 1.72^m, O<0<1 

(l+m) w -"-|A n | = 1.72 07im + l-- |A*|, 

and therefore by (82) (l + m)"-*— \A n \^1.72dnm+m/M. 

For n^4 we have in the case (82), since the inequality (81) is not satisfied, 

±<JL 
M^l.5 

and, therefore, finally 

(l+m) n -"-\A n \^2.39nm (*>=1, ...,n) (83) 

On the other hand, since Af<1.5/n, (l+2W)'- 1 ^(l + 1.5/7i)*<e 1 - 5 =4.481689, and, therefore, 
from (72) and (73) 

t-tr ^ 10.72n |£| („=1, ...,ti; |A.| J^l-Jj 7i>3^, 

[A- 1 -^- 1 !^ 10.7271t£| (p=l,»; |A w |^l-|p ti>3^ (84) 

To obtain a lower bound, we take in (74) v=l. We obtain, since for ti>3, 0<|A n | £ 1, 

K»| Kl ~~ l A «l 
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and, therefore, 



m 



lA- 1 -^ ' ~%^ 0.75u j—r &>=!.« ;n£4). 



(85) 



To obtain the exact value of [A" 1 — A ( £ ) ~ 1 \ p if (75) is not satisfied, we return to the inequali- 
ties (78) which are equivalent to (77). 

Denote by n the smallest integer between 1 and n, such that we have 



r» n -i ^ |Aj>o- w 



(lSn»S»)- 



(86) 



The parts of this inequality implying a or <r n must be disregarded, that is to say, this inequality 
reduces to |A w |>o-! for n =l and to er„_i^|A n | for n =n. Then we see at once that 
max (t v —tf)=t nQ —t { o } , and therefore 



— / — / (0) 



(1 +M) n o~ 1 (1 +m) n - n o-\A n 



\&*\ 



(p=l,oo). 



For m=M we have %=1, and therefore by (23) 

,(o)-ii _ (ft-l)m 



AT l -A 



(m = M<^ I ) 



(87) 



(87°) 



1— (n — l)m 

7. Bounds of A-'-A^- 1 

In order to obtain the theorem corresponding to E for a determinant A satisfying condi- 
tions (56), (57) of theorem B, we prove first the following important lemma, which generalizes 
considerably the relation (52°) of lemma I. 

Lemma III. Consider n positive numbers a i} . . ., a n and 2n 2 —2n nonnegative numbers 
flV, e M „ (n^v, tx,v = l } . . ., n) , such that the matrices 



M* 



«1 — ^12 — €12 

— m 2 \ — € 2 l «2 



k. — *W»1 6 «1 m n2 — C«2 

^ «l — ^12 

-m 2 i a 2 

M (0) = 



Win € lw 

^2w — €2 W 



«n J> 



(88) 



— ^2« 



k.— w„i — m n2 
are M-matrices. 

Consider n constants A^ (p—l, , . ., n), such that 

lA^l^a, (m=1, . . ., n) 



(89) 



(90) 
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and 2ri 2 —2ri constants a ilV) b^ (v^v; p,? 5 !, . . ., n) such that 

and form the two matrices 

r A x bi2 ... 61 ^ 



(91) 



A= 



621 ^2 



^4i a 12 

#21 ^±2 



A m = 



^flnl Unl 



b ir , 



(92) 



AnJ, 

a 2 



j&nj* 



(93) 



Then we have 



A- l -A i0) - l <^M- l -M^- 1 . 



(94) 



Proof: We write M=P—T, M (0) =P— T (0) , where the matrices T, T (0) have zeros 
in the main diagonal and off the main diagonal, respectively, the nonnegative elements m^+^j 
ra M „ and where P is the diagonal matrix 



Oil 



a 2 











<w 



We can develop the inverse of P— T=P(E—P l T) (E is the unit matrix) in the following way: 

(P-T)~ l = i2{P- l Typ-\ (96) 

K = 

The convergence of this development and the validity of (96) follow easily from the fact 
that the determinant of the matrix P—tT does not vanish for \t\ ^ 1, which follows immediately 
from lemma I. The corresponding development holds also for (P— T (0) ) _1 , and we obtain 
therefore 



M- l -M^- l ^^[{P- l ry-{p- l T^y]P-K 

K = l 



(97) 
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The elements of P T are here (m Mi , +«„,)/«„ or zeros, and those of P 1 T° are m,,„/a„ or zeros; 
therefore all elements of the matrices 



(p-i2y_(p-i r <o>)« ( K =l,2, . . .) 



(98) 



are polynomials in ra M „/a M and e^/c^ with nonnegative coefficients. Denote by (? the diagonal 
matrix 



(A, 



Q- 











A n> 



nd write ^L=#— S, A {0) = Q— # (0) , where the elements of S are b^ or zeros, and those of 
S (0) , a?, or zeros. Since, therefore, we have 



we have 



Q-'S^P-'T, Q-'S^^P-'T™, 



K = K = 



A- 1 -A™- 1 =J}[(Q- 1 S)'--(Q- l S {0) y]Q- 1 > 



(99) 



But now the elements of (Q~ 1 S) K —(Q 1 S {0) ) K are obtained from those of (98) in substituting 
there a^/A^ instead of m^/a^ and (b^—a^/A^ instead of e M „/a M . Further, we have by (90) 
and (91) 



.m u 



O fiV O/j^p 



and therefore, since the coefficients in (98) are as already mentioned not negative, 

{Q- l sy~-(Q- l s {0) y^:{P- l TY-{P- l T 0) Y 



and 



^[{Q- l S) K -{Q- l S^Y]<^Z[(P- l TY-{P- l T^Y]. 



(100) 



But from (100) by virtue of (90) it follows that the development (99) is majorized by (97) 
and our lemma is proved. 

Under the conditions of theorem B in section 4, if the inequalities (56) and (57) are satisfied, 
we can apply lemma III in replacing the matrix M by A n , M (0) by A®\ We obtain 

F. Under the conditions of theorem B, if the inequalities (56) and (57) are satisfied, we have 



| 4 _ 1 _ A( 0)-i|^| A -i_ A (o)-i| p=max(1+M) 



^ (1+m)"-"-)^! 



(P=1,co), (101) 



where the values and the estimates for the right-hand expression in (101) are obtained from the 
formulas (79) to (87°). 

Lemma III can be applied to many problems similar to that solved by theorem F. For 
instance, in the theory of the solutions of linear equations, the following problem has to be 
dealt with, although its complete discussion is usually avoided: 
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Consider a ' 'triangular" system of linear equations 

• 



(102) 



where the coefficients a M „ are only approximate values to the "true" values b^. Suppose that 
we have generally 

K-a^l^e (jr^ M ); (103) 



how far is the solution of the system 



7 I V ILV<*'V Up 



(104) 



influenced if the b^ are replaced by a^l 

If we denote the matrix of (102) by A i0) and that of (104) by A, the question can be 
answered by giving estimates of \A~ l —A m ~ l \ p (p=l, °°). Suppose that we have generally 



ow=i, W^Afji (m>^; m=i, • • • ,w— 1) 
and consider the matrix 



(105) 



A(Mj, 



^.-, 



1 



-M 







1 



(106) 



where all elements to the left of the main diagonal are zeros, and all elements to the right of 
the main diagonal in the /ith row are equal to — M M . 

We obtain from lemma III at once the majorization 

^_ 1 _^ (0) _ 1<A(Mi + e? ? M w _ 1 + e)- 1 -A(M b . . ., M n _ x )-\ 

To obtain the inverse of the matrix A(Mi, . . ., M n _i), consider the system of linear equations 
x ll —M ft (x fi + 1 + . . . +x n )=u ll (m=1, • • ., n). (107) 

To solve it explicitly, we put 

* M =ZX (m=1, . • ., »), (108) 

1 + M^iV,, AyV M+1 . . . N n =P, ( M =l, . . ., n-1) 2VV=1, P„=l. (109) 

Then (107) becomes 

Xv—M^+i+Ur (m=1, . . ., ft), (U0) 

s^A^h-i+^m (m=1, . . ■., n), (111) 

where s n+1 =0. Dividing (111) by P M , we obtain 
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and therefore 



introducing this in (110) we have 



^> u K . 



5 M *V-jL-I p 7 

k=/j. * k 



n m 

K=fl-\-l ±K 



(112) 



We obtain therefore for the inverse of our A (Mi, . . ., M„_i), if E denotes the unit matrix, 

A(M 1; . . ., M n ^=E+DT, (113) 



where D is the diagonal matrix 



D-- 



M,P 2 



M 2 P 3 







and 7" the triangular matrix 



T= 



The expressions 



J_ J_ 
F 2 P 3 



° P, 











M„_,P M 



Pn 

\_ 
Pn 

\_ 
Pn 



J. 



A(A/ 1 + e, . . ., M K _ 1 + e)-'-A(M 1 , . . ., M,^)" 1 



(114) 



(115) 



(116) 



obtained from (113), (114), (115) are rather unwieldy; however, we shall obtain for its norm 
corresponding to p = l (cf. 66, a) a very simple and elegant expression. 

Indeed, if we take the sum of the elements in the *>th column of A(M if . . ., M w _i) -1 and 
denote it by t' v , we have 

If we now write out t' v+1 , we obtain 

Ci=l+rS^ n 1 (1+M K )1(1+M V )+M V ; 
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and in comparing this with the expression for t\ we see that we have C+i s -K&QiJ and 
therefore 

t' v =n(i+M k ). (117) 

We obtain now for the sum of the elements in the *>th column of (116) 

t v =n\i+M K +e)-n\i+M K ) (F=l,..-. f ii). 

In multiplying this by l+M„+c, we obtain 

n(i+M K +e)- n(i+M,)-en(i+^), 

and comparing this with t v+u 

t v+1 = (l+M v + e)t v + eU(l+M K ). (118) 

*c=l 

Therefore, /„ is monotonically increasing with v, and we see that the norm of (116) cor- 
responding to p = l has the value 

n\l+M,+e)- n n\l+M 9 ) (119) 

and obtain therefore 

\A- 1 -A^- l \ 1 ^n\l + M v + e)- n Tl\l+M v ) (120) 

p=i p=i 

as the solution of our problem. In applications it may be better to use the recurrence formula 
(118). If all M M have the same value M, the expression (113) coincides with that obtained in 
section 5 for A" 1 for m = 0. But in this case we see from (48) and (49) that the row sums run 
through the same values as the sums in the columns. We obtain, therefore, in this case the 
expression 

(M+l + c) w - 1 -(M+l) w - 1 =(^-l)€(M+l4-<9€) ? O<0<1, 

as the norm of (116), both for p= 1 and p= c° . 

8. Linear Systems With a Nearly Triangular Matrix 

The results of the preceding sections give the means to discuss the following problem 
concerning the system (1) in the introduction under the conditions (2) and the "triangular" 
system 

n 

X>mv^=2/m (»=!,.. .,n), (121) 

with the matrix A m . In discussing this problem we can obviously assume that 
^mm^I (m=1j . . . ,n). Then the difference between the solution of (1) and that of (121) is 
given by the vector (A~ 1 —A (0) ~ 1 )t] j rj=(y h . . . ,# n ), and the norm of this vector corresponding 
to one of the indices p=l, «> does not exceed |^i -1 — ^1 (0) ~ 1 | 7 ,It7| / , (p = l,co) f and can indeed for 
suitable choice of the vector tj attain this limit. Therefore, the norms \A~ X — A {0) ~ l \ p measure 
the error committed in replacing the system (1) by (121). 

For an e>0,M>0 being given, how small must m^>0 be taken in order that we have 

\A- l -A*>-%£* 0=l,oo)? (122) 
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If we introduce the quantities |A ra | and 5 corresponding by (22), (23), and (27) to m and 
M, we obtain from (79) and (94) the condition 



(1+M) ,-i_A_ =€ 



<5< 



(l+M^ + e 



(w*i-5) 



(123) 



as long as the condition (75) is satisfied and therefore certainly as long as M = 1.5/n (n = 4:) 
On the other hand, we have by (30) 

8 



m-- 



M " + e M 



and from (123) and (124) 



M n (l+.M)»- 1 +(M n +-^)« 
We see that it will be sufficient to take 



M n {l + M)^+(M H +-^y 



(ki*i-£> 



(wi"-5) 



(124) 



(125) 



(126) 



to solve our problem, for instance, if we have M^l.5/n (n^4) or M^O.5321 (n,=3). For 
small values of e obviously, only the first term in the denominator is essential, and we have 



m ~K(nM); K ( n > M )= Mn (l+M) n - 1 ' 



(127) 



Tables 1 and 2 give the values of K(n,M) for a set of integer M from 1 to 10 and some values 
of M^>1.5/n. We have obviously 

K(n } M)<j~^ 2 ' (128) 

Table 1. K(n,M) 



M 


n = 3 


71=4: 


n=5 


n = 6 


n = 7 


n = S 


71 = 9 


n=10 


1 
2 
3 

4 
5 

6 
7 
8 
9 
10 


.0625 
.0111 

. (2)347 
. (2)143 
. (3)694 

. (3)378 
. (3)223 
. (3)140 
. (4)926 
. (4)636 


.0114 

. (2)103 

. (3)193 
. (4)526 
. (4)182 

. (5)736 
. (5)336 

. (5)168 
. (6)903 
. (6)515 


. (2)240 
. (3)106 
. (4)116 
. (5)206 
. (6)498 

. (6)149 

. (7)522 
. (7)207 
. (8)900 
. (8)424 


. (3)548 
. (4)115 
. (6)719 

. (7)821 
. (7)138 

. (8)304 
. (9)815 
. (9)255 
. (10)900 
. (10)351 


. (3)130 

. (5)126 
. (7)448 
. (8)328 
. (9)383 

. (10)619 
. (10)127 
. (11)315 
. (12)900 

. (12)290 


. (4)316 
. (6)140 
. (8)280 
. (9)131 
. (10)106 

. (11)126 
. (12)199 
. (13)389 
. (14)900 

. (14)239 


. (5)778 
. (7)155 
. (9)175 
. (11)524 

. (12)295 

. (13)258 
. (14)311 
. (15H80 
. (16)900 
. (16)198 


. (5)193 
. (8)172 
. (10)109 
. (12)210 
. (14)821 

. (15)526 
. (16)486 
. (17)592 
. (18)900 
. (18)164 
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Table 2. K (n, M) 



M 


n = S 


n = 4 


ti = 5 


71 = 6 


n = 7 


n = 8 


7i = 9 


71=10 


. 1 
, 2 
.3 
.4 
. 5 

. 6 
.7 
.8 
. 9 
1.0 

. 1 
. 2 
.3 
. 4 
.5 

. 6 

. 7 

.8 

.9 

1.0 


~254~~ 

. 181 
. 134 
. 102 
.0789 
.U625 


"YlY" 
.0718 

.0464 
. 0313 
.0218 
.0156 
.0114 


.~0866~ 

.0438 

.0241 

.0141 
. 00864 
. 00548 
. 00359 
. 00240 


.0399"" 
.0180 

. (2)891 

. (2)470 
. (2)260 
. (2)150 
. (3)894 
. (3)548 


.0196"" 

. (2)788 
. (2)349 

. (2)165 
. (3)825 
. (3)431 
. (3)233 
. (3)130 


.0828 " 

. 0101 

. (2)359 
. (2)142 

. (3)602 
. (3)270 
. (3)127 
. (4)623 
. (4)316 


. 0197" " 
. (2)533 
. (2)169 
. (3)592 

. (3)224 
. (4)902 
. (4)382 
. (4)169 
. (5)778 


"0Y2T" 

. (2)289 
. (3)809 

. (3)252 

. (4)848 
. (4)305 
. (4)116 
. (5)462 
. (5)193 


7i=15 


n = 20 


tt = 25 


7i = 30 


7i=35 


7i = 40 


n = 45 


7i = 50 


.0157 
. (2)137 
. (3)167 
. (4)242 
. (5)399 

. (6)729 
. (6)146 
. (7)317 
. (8)743 
. (8)186 


. (2)439 
. (3)188 
. (4)112 
. (6)809 
. (7)681 

. (8)658 
. (9)721 
. (10)886 
. (10)121 
. (11)182 


. (2)138 
. (4)281 
. (6)793 
. (7)277 
. (8)118 

. (10)597 
. (11)357 
. (12)248 
. (13)197 
. (14)178 


. (3)469 
. (5)439 
. (7)570 
. (9)957 
. (10)204 

. (12)543 
. (13)177 
. (15)695 
. (16)322 
. (17)173 


. (3)166 
. (6)697 
. (8)413 
. (10)331 
. (12)354 

. (14)494 
. (16)879 
. (17)195 
. (19)525 
. (20)169 


. (4)604 
. (6)112 
. (9)299 
. (11)114 
. (14)613 

. (16)449 
. (18)436 
. (20)545 
. (22)856 
. (23)165 


-. (4)224 
. (7)180 
. (10)217 
. (13)395 
. (15)106 

. (18)409 
. (20)216 
. (22)153 
. (24)140 
. (26)162 


. (5)841 
. (8)290 
. (11)157 
. (14)137 
. (17)184 

. (20)372 
. (22)107 
. (25)428 
. (27)228 
. (29)158 



The bound in (126) is obviously the "best" under the condition (75), save that the factor 
1/M of e in the denominator could be replaced by an (unknown) fraction of it. 

If M<^1.5/n, the use of the general formula (87) is very cumbersome. We can, however, 
obtain a good working limit for m in the following way. If the inequalities (2) are valid for 
an M<C.l.5/n=M (n) , then they are also satisfied if M is replaced by M in) , but then the limit 
for m obtained for M {n) is also sufficient for our M. We obtain therefore in this case the 
sufficient condition for (122) in the form 



ra<ran=- 



Mi n) (l+M (n) ) n - 1 + 



( M "+£> 



(129) 
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Table 3 gives the values of K(n,M {n) ) for n=l,2, . . ,, 50. 

Table 3. M = M^ = 1.5/n 



n 


(l + M)"-nM-l 


K(n,M) 


n 


(l + A0 n -nil£-l 


K(n,M) 


3 


.871 


. 255 


27 


1. 805 


. 00755 


4 


1. 074 


. 134 


28 


1. 811 


. 00723 


5 


1. 213 


. 0866 


29 


1. 817 


. 00694 


6 


1.315 


. 0623 


30 


1.822 


. 00667 


7 


1. 393 


. 0480 


31 


1. 827 


. 00642 


8 


1.454 


. 0387 


32 


1. 831 


. 00619 


9 


1.504 


. 0323 


33 


1.836 


. 00597 


10 


1. 546 


. 0276 


34 


1. 840 


. 00577 


11 


1. 580 


. 0240 


35 


1.844 


. 00558 


12 


1.610 


. 0213 


36 


1. 848 


. 00540 


13 


1.635 


. 0190 


37 


1. 851 


. 00524 


14 


1. 658 


. 0172 


38 


1. 854 


. 00508 


15 


1. 677 


. 0157 


39 


1. 857 


. 00493 


16 


1. 695 


. 0144 


40 


1. 860 


. 00480 


17 


1. 710 


. 0133 


41 


1. 863 


. 00466 


18 


1.724 


. 0124 


42 


1.866 


. 00454 


19 


1. 736 


. 0116 


43 


1. 869 


. 00442 


20 


1. 748 


. 0109 


44 


1. 871 


. 00431 


21 


1.758 


.0102 


45 


1. 873 


. 00420 


22 


1.768 


. 00965 


46 


1. 876 


. 00410 


23 


1. 776 


. 00914 


47 


1. 878 


. 00401 


24 


1. 784 


. 00869 


48 


1. 880 


. 00391 


25 


1.792 


. 00827 


49 


1.882 


. 00383 


26 


1. 799 


. 00789 


50 


1. 884 


. 00374 



The expression for m' can be written in introducing the value of M ( n n) as 

^'o+¥no+¥)"-")+(o+¥)"-'-«y 

We observe now that for any positive a and positive x the expression (1 + x/n) n ~ a mono- 
tonically increases and tends to e x if the positive n increases monotonically to oo . Indeed, if 
we put u=l/n, take the logarithm of this expression, differentiate it with respect to u, and 
multiply by u 2 , we obtain 



x(u — au 2 ) , 7- -. x 
— ~ -—log (1 +ux) 



but this expression vanishes for ^=0 and decreases for positive u, since its derivative is 

— xu 
(H 

We see, therefore, that 



(1+ux)' 



(2a-\-x-\-axu). 



o +¥)">'■•. o+^y^- 



and our bound for m can therefore be replaced by 

3 € ^ € 

2n ' e U5 (e U5 - 2. 5) + (e 1 - 5 - 1.5)6^(5.93 + 1.99eV 
we can replace therefore the condition m^m' by the simpler condition 

€ 1 



m^m = 



n 6 + 2e 



(130) 
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On the other hand, we obtain at once the solution of our problem in the case m = M< l/(n— 1) 
from (87°). 

m^m ^^--^— (m = M<~ ) —\ (131) 

This is a "best" condition, whereas the condition (130) can still be improved. We see 
however, in comparing (130) with (131), that the bound in (130) cannot be improved by a 
factor greater than 6. 

From the condition (130) we can finally derive the inequality (7) of the introduction. 
Indeed, if positive m <ll/2n is given, we can solve (130) with respect to e, 

6nm 



1 — 2nm 
and obtain therefore, in applying (130) for m=m , the inequality 

\A-*-A*>-*L< ■ 6nm ° , 
1 {p=z l-2nm 

which holds for all positive ra <l/2?i and gives the inequality (7) if we replace here m by m. 
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